Cosmological Information from Lensed CMB Power Spectra 
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Gravitational lensing distorts the cosmic microwave background (CMB) temperature and po- 
larization fields and encodes valuable information on distances and growth rates at intermediate 
redshifts into the lensed power spectra. The non-Gaussian bandpower covariance induced by the 
lenses is negligible to I = 2000 for all but the B polarization field where it increases the net variance 
by up to a factor of 10 and favors an observing strategy with 3 times more area than if it were Gaus- 
sian. To quantify the cosmological information, we introduce two lensing observables, characterizing 
nearly all of the information, which simplify the study of non-Gaussian impact, parameter degen- 
eracies, dark energy models, and complementarity with other cosmological probes. Information on 
the intermediate redshift parameters rapidly becomes limited by constraints on the cold dark matter 
density and initial amplitude of fluctuations as observations improve. Extraction of this information 
requires deep polarization measurements on only 5-10% of the sky, and can improve Planck lensing 
constraints by a factor of ~ 2 — 3 on any one of the parameters {wo,Wa, ^k, Y1 ™-^} with the others 
fixed. Sensitivity to the curvature and neutrino mass are the highest due to the high redshift weight 
of CMB lensing but degeneracies between the parameters must be broken externally. 
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I. INTRODUCTION 

Primary cosmic microwave background (CMB) 
anisotropy from recombination has proven itself to be 
a veritable gold mine of cosmological information. One 
of the most important secondary signals that should be 
detected by upcoming cosmic microwave background 
experiments is the distortion to the temperature and 
polarization fields due to gravitational lensing by the 
large-scale structure of the universe (see 3 for a recent 
review). Lensing distortions add cosmological informa- 
tion on parameters such as curvature, neutrino masses 
and dark energy that change the expansion and growth 
rate at intermediate redshifts {z <5). 

This distortion in real space couples power in har- 
monic space and hence introduces non-Gaussianity into 
the CMB temperature and polarization fields. Beyond 
power spectra, this non-Gaussianity is a source of infor- 
mation in that it allows direct reconstruction of the con- 
vergence field 0,11,0, IE Hi' On the other hand, for pur- 
poses of extracting cosmological information from lensed 
power spectra as considered here, this non-Gaussianity is 
largely an impediment as it makes power spectrum esti- 
mates covary across a wide range of multipoles. 

The purpose of this paper is twofold. First, we calcu- 
late the full non-Gaussian covariance between all combi- 
nations of temperature and polarization bandpowers in 
the lensed CMB. This extends previous work in which 
the temperature 0, 13 and _B-mode polarization covari- 
ance Q were calculated separately. Second, we present 
a general framework for studying the extra information 
on cosmological parameters that lensed CMB spectra 
supply, with particular attention to the impact of non- 
Gaussianity. 

Previous works have noted that the lensed CMB signal 
may be used to study the dark energy 0, 0, ^| and 



neutrino mass [T^l- These studies did not compute the 
non-Gaussian covariance but assumed either that the in- 
formation is encoded in the unlensed primary CMB and a 
reconstruction of the lenses or by approximating the non- 
Gaussian covariances with a degradation factor from . 
Our results lend support to these analyses. We also study 
the sensitivity of lensing to curvature and find that future 
CMB measurements can provide interesting constraints 
on it. 

This paper is organized as follows. In we compute 
non-Gaussian contributions to the covariance between 
all lensed CMB temperature and polarization bandpow- 
ers. We then describe in illlll how this non-Gaussian co- 
variance propagates into Fisher matrix parameter fore- 
casts and present formal bounds on its impact. In illVI 
we define two parameter independent observables which 
contain essentially all information from the lensed CMB 
and discuss their relationship to distance and growth as 
well as their degeneracy with parameters that control the 
matter power spectrum. Armed with this general frame- 
work, we show how constraints on these observables can 
be interpreted in the context of common parameteriza- 
tions of the dark energy and dark matter in Wl Finally in 
ilVIl we show how future CMB surveys can be optimized 
for sensitivity to the lensing observables. We conclude in 
j|VIII and briefly address the issues of goodness-of-fit in 
Appendix ^ and scaling with the fiducial cosmology in 
Appendix ^ 



II. LENSED POWER SPECTRUM 
COVARIANCE 

In this section, we compute the non-Gaussian covari- 
ance between all CMB temperature and polarization 
bandpowers to lowest order in the lensing power spec- 
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trum Cf'^. The results (Eqs. (O, |THll, and ^) will 
be foundational in subsequent sections, as they will per- 
mit the effects of non-Gaussianity to incorporated into 
parameter forecasts. However, the details of the calcula- 
tion will not be needed, so the reader may wish to skip 
this section on a first reading. 

First, we recall some preliminaries concerning lensed 
CMB fields. We work in the flat sky approximation; 
we will see f ^III B|l that non-Gaussian covariance only 
becomes important when combining bandpowers over a 
wide range of ^, so that all-sky corrections from the dis- 
crete nature of I should be negligible. The lensed CMB 
temperature T(x) and unlensed temperature r(x) are 
related by 



T(x) =T(x + V</.(x)) 



(1) 



The projected potential is given by the line-of-sight 
integral: 

where D — J dz/H is the comoving distance along the 
line of sight, Dg denotes the comoving distance to the 
surface of last scattering, and 



DAiD) 



1 



sinh JQkH^D 



(3) 



is the comoving angular diameter distance. 

Polarization fields are lensed in the same way; the 
lensed Stokes parameters Q(x), C/(x) and unlensed ver- 
sions (5(x), C/(x) are related by 



g(x) = g(x + v0(x)), 

[/(x) = J7(x + V0(x)). 
The Fourier versions of Eqs. QJ, l@J are 



(4) 



T(i) = r(i) - 

Eil) = ^(1) 



(2^)2 
dH' 
{2n) 



WTil,l')T{V)(f>{l-V) + 0{^^), 
^WE{l,l')E{l')cf,{l-l')+Oi<l>'), 



B{\) 



21/ 



;WB{ll')E{l')cf>{l~l')+0{ 



(27r)2 

where the kernels are defined by 



(5) 



WtM - -[l'-(l-l')], 

We{1,1') = -[1' -(1-1')] cos 2(^1 -(^lO, 

Wb{1,1') = [1' -(1-1')] sin 2(^1 -(^lO, (6) 

and the E and B fields are related to the Stokes param- 
eters as 

[E{1) ± tB{l)] - [g(l) ± tU{l)] exp(T2z^i) . (7) 



Here (p\ is the angle between 1 and x axis. We have 
assumed that the unlensed_B = such that the observed 
B field is generated from E by lensing alone [THj l. 
We define ideal, noise-free bandpower estimators by 

W = dH(J^\x*{\)Y{\)^ (8) 



27r 



where XY e {TT, EE, TE, BB}, A is the survey area in 
steradians, and 



d'l 



(9) 



lei 



is the Z-space area of band i. We define power spectra as 
usual 



such that 



Ar = {AD 



A6i,vCi 



d'l 



XY 



2tt 



XY 



(10) 



(11) 



(As a technical point, when propagating C[ to lensed 
power spectra, we use the all-sky correlation function ap- 
proach of 0, for consistency with CAMB.) 

Now let us consider the covariance of these estimators 
in non-overlapping I bands. We split the bandpower co- 
variance into Gaussian and non-Gaussian pieces: 



-,XY,ZW 



XY.ZW 
ij 



(12) 



Here and below we will use the short hand notation / to 
denote a unique bandpower specified by the Z-band i and 
the power spectrum XY. The Gaussian piece is given 
by: 



^XY.ZW 



Aa? 



dH fil 



(13) 



where the power spectra which appear are lensed. In the 
presence of instrumental noise, the power spectra in this 
formula are replaced as 



XY 



XY 



(14) 



where the noise power spectra for white detector noise 
with a Gaussian beam are given by 



XX 



iln2 



\ TqMB ) 



(15) 



for XX G TT, EE, BB and vanishing for other spectra. 
We will also take Ag = As = Ap. 
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In 13], we computed the non-Gaussian piece A/^ 
for the case of two BB bandpowers: 



XY,ZW 
ij 



BB.BB 2 f 2, f ,2i f '^^^ 



(16) 



where 



= wuh,muh^){crfctf^,ct;_,, (i?) 

with Cf^^ ^= as the unlensed power spectra. 

Here, we consider two additional cases. First, we com- 
pute the non-Gaussian covariance of one BB bandpower 
with a bandpower Af ^, where X,Y e {T, E}: 



■XY.BB 



Aa.aj Ji^^i " Ji^gj ^ (27r)2 ' 

xiyB(i„i,)'Qf^Qf^Qt-i, • (18) 



Second, we compute the covariance of two bandpowers 
AP, Af^, where X,Y,Z,W e {T,E}: 



^XY,ZW _ 



/ XY.ZW , ZW.XY , qXY.ZW 



XY. 



, oXY.WZ , pYX,ZW , aYX,WZ\ 

+(3uu +^1.1, ) 



Here, 



«M/^^= W^X(1., h)Wy{h, l,)(Qf ^Q^^) , 

AY'^^^ Wxih, WWz{h,h)C,^^Cl^ . (20) 



paper, we use a fiducial model consistent with the third- 
year WMAP lig data: 

{nbh'^,nch'^,T,\n S(;,n„r} = {0.0223,0.104,0.088, 

-10,0.951,0}, (23) 

{r!DE,^^i^,^m„u;o,u;a} = {0.76, 0, 0.061eV, -1, 0} . 

The first set represents parameter that control the intrin- 
sic power spectra from recombination whereas the second 
set represents the intermediate redshift parameters of in- 
terest to lensing. Here 5,^ oc vG4 is the amplitude of 
initial curvature fluctuations at fc = 0.05 Mpc~^ and r is 
the tensor to scalar ratio. The dark energy (DE) equation 
of state is parameterized as 



w{a) = Wo + (1 — a)wa ■ 



(24) 



When Wa = we will use the variables w and wq inter- 
changeably. 

The most important part of the non-Gaussian covari- 
ance is between two BB bandpowers, where the ana- 
lytic and Monte Carlo results agree well. The second 
most important contribution is the correlations between 
EE and BB which appear at the 10-30% level. Here 
the agreement between Eq. ()18|l and the Monte Carlo 
results for the correlation are at the 10-20% level indi- 
cating that higher order contributions are not entirely 
negligible. However they represent a small correction to 
a small correlation and we neglect it throughout. In the 
remaining parts of the covariance matrix, including the 
entire {TT, TE, EE} covariance, the non-Gaussian con- 
tributions are small. 

Since the non-Gaussianity manifests itself as a small 
correlation across a wide range of multipoles, it is only vi- 
sually apparent when combining the multipoles into large 
bands as in Tab. In practice, in the following sections 
we take bins of A^ = 10 out to Ij^ax = 2000, and all 
results in the paper are robust to binning more finely. 



III. PARAMETER FORECAST FORMALISM 



Taken together, Eqs. (|16|l . (|18|l . H19|l constitute a com- 
plete calculation of the (4iVband)-by-(4iVband) covariance 
matrix. 

In Table m we show the correlations 



R 



IJ 



and variance degradation factors 



D' 



(CO) 



II 



(21) 



(22) 



for aU combinations of {TT, TE, EE, BB} in two large 
I bands, with Monte Carlo results from 10^ simulations 
shown for comparison. Here represents the Gaussian 
contribution to the covariance matrix. Throughout this 



In this section, we review the Fisher matrix formalism 
for forecasting parameter errors in the presence of non- 
Gaussian errors in illll Al We then apply this formalism 
to place formal bounds on the impact of non-Gaussianity 
in TObI 



A. Fisher matrix 

The Fisher matrix provides a useful way of assessing 
the impact of the non-Gaussian bandpower covariance 
on parameter estimation. Even for cases where the like- 
lihood function cannot be evaluated directly, it can be 
approximated as the linear propagation of errors from 
bandpower space to another parameter space Pa ■ In Ol 
we gave a complete calculation, to lowest order in Cf'^, of 
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TTio 


TTm 


TEio 


TEhi EE\o EEhi 


BBio 


BBm 


D 




1 


0.007 


-0.053 


0.001 0.074 0.001 


0.025 


0.009 


1.007 (1.012) 


TThi 


(0.008) 


1 


0.001 


-0.312 0.003 0.089 


0.014 


0.025 


1.020 (1.019) 




(-0.055) 


(0.002) 


1 


0.003 -0.098 0.001 


-0.036 


-0.010 


1.000 (1.000) 


rShi 


(0.001) 


(-0.311) 


(0.004) 


1 0.001 -0.306 


-0.049 


-0.086 


1.010 (1.011) 


EEio 


(0.076) 


(-0.001) (-0.096) 


(0.001) 1 0.004 


0.316 


0.137 


1.012 (1.011) 




(0.002) 


(0.090) 


(0.002) 


(-0.311) (0.003) 1 


0.137 


0.283 


1.039 (1.039) 


BBla 


(0.022) 


(0.027) 


(-0.048) (-0.030) (0.311) (0.117) 


1 


0.754 


4.323 (4.416) 




(0.005) 


(0.039) 


(-0.021) (-0.067) (0.132) (0.262) (0.754) 


1 


7.595 (7.619) 



TABLE I: Bandpower correlations and variance degradation factors D (see Eq. 1221 ') using two bands in /: (lo)= [100, 1000] and 
(hi)= [1000,2000]. Upper diagonal values were calculated analytically using the lowest-order expressions Eqs. I16i . 1181 . 1191 : lower 
diagonal values (in parentheses) were calculated from 10^ Monte Carlo simulations. The large correlations between TE and {TT, EE} 
are dominated by the Gaussian contribution. 



the (4iVband)-by-(4iVband) covariance C^"^ between bands 
/ and J specified by the power spectrum combination and 
the / range (see Eq. (|12|l ). In terms of this covariance, 
we define an approximate Fisher matrix as 



i^a/3 = ^(a„A^)(c-i),j(a;3A^), 



(25) 



IJ 



where a, /3 run over a basis set of directions in parameter 
space. In this section we will use upper indices to denote 
quantities that transform as a contravariant tensor under 
a re-parameterization and lower indices for those that 
transform as a covariant tensor. For example, Eq. H25(l 
represents the transformation of the inverse covariance 
matrix from the bandpower space to the parameter basis 
space. 

Given this transformation, the inverse of the Fisher 
matrix can be thought of as an estimate of the covari- 
ance matrix of the basis parameters. As such it gives the 
variance of any hnear combination of basis parameters as 



Var(7r) =^(5„^)(F-i)"'3(a^7r) 



(26) 



Q/3 



As a special case, if tt corresponds to a basis direction a, 
then the marginahzed uncertainty is given by the diago- 
nal element (F^i)"". 

The Fisher matrix quantifies the local curvature of 
the likelihood function in the parameter space. Fisher 
forecasts therefore suffer from several problems (see e.g., 
01 El) J especially in the presence of nearly degenerate 
parameter directions. If the derivatives in Eq. (|25|l are 
not constant across the extent of the degeneracy, this 
curvature is also not constant and hence confidence re- 
gions in parameter space are not well approximated by 
ellipsoids in either their shape or extent. Fisher matrix 
forecasts should only be interpreted as confidence lim- 
its on parameters for well-constrained directions in the 
parameter space and as a tool to expose parameter de- 
generacies. Finally even if the parameter derivatives vary 
significantly only outside of the predicted error ellipsoid. 



Fisher matrix forecasts still depend on the choice of the 
fiducial model. These points should be kept in mind when 
considering the parameter forecasts in the next two sec- 
tions. 



B. Formal Bounds on Non- Gaussian Impact 

Before turning to parameter forecasts in specific pa- 
rameter spaces, it is instructive to quantify general 
bounds on the impact of non-Gaussianity. One of the 
main results of this paper is that non-Gaussian power 
spectrum covariance is essentially negligible when con- 
sidering lensed {T, E} alone at Zmax = 2000 (c.f. for 
higher Zmax)- Beyond this Imax other secondary sources 
of temperature and polarization will likely prohibit the 
extraction of cosmological information. 

To state this in a precise way, we introduce Karhunen- 
Loeve (KL) eigenvalues between the non-Gaussian and 
Gaussian bandpower covariances. The KL eigenvalues 
Xk and eigenvectors are defined by 



A 



(27) 



where C, C*^ denote the bandpower covariances with and 
without non-Gaussian contributions. With cosmic vari- 
ance limited {TT, TE, EE} bandpowers to Z^ax = 2000, 
we find that each Xk is between Amin — 0.94 and A„iax — 
1.08 in the fiducial cosmology. The exact values depend 
on the normalization of the power spectrum but remain 
around unity for all reasonable variations. 

These eigenvalues limit the excess variance on parame- 
ter errors from non-Gaussianity. More precisely, we now 
prove that for any cosmological parameter tt, the ratio of 
uncertainties with and without including non-Gaussian 
covariance satisfies: 



An 



Var(^ . , 
- Var°(7r) " 



(28) 



This inequality holds even after marginalizing any set of 
additional parameters. 
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The first step in the proof of Eq. (|28|l is to note that in 
a basis consisting of KL eigenvectors in the bandpowers, 
both the non-Gaussian and Gaussian bandpower covari- 
ance matrices are diagonah 



(29) 



(More generaUy, this holds true for symmetric matrices 
as long as one of them is positive-definite.) Now consider 
any estimator E which is hnear in the bandpowers. In 
the KL basis, it can be written 



(30) 



K 



Combining Eqs. (|29|l and H3()() . the ratio between the 
estimator variance calculated with and without non- 
Gaussian covariance satisfies: 



Var(£) _ Y.k>^k{£k? 



Var^(£) Y.k{^k) 
In this form, it is seen that 

^ Var(g) 
Var^(£) 



(31) 



(32) 



Next we observe that in the Fisher approximation, the 
estimator for any cosmological parameter tt depends 
linearly on the estimated bandpowers. (This is still 
true if additional parameters are marginalized, although 
marginalization will change the optimal estimator.) We 
denote the estimator which is optimal with the non- 
Gaussian covariance included by tt and that without ttq . 
Then: 



A 



< 



Var(7f) ^ Var(7f) 



min ^ c , 

Var°(^) 



Var°(^G) 



< 



Var(7fG) 
Var°(^G) 



< An 



(33) 

where we have combined Eq. H32|l with the inequalities 
Var°(7fG) < Var°(7f) and Var(7f) < Var(7fG), which fol- 
low from the optimality of each estimator. This com- 
pletes the proof of Eq. . 

The KL eigenvectors also illuminate the nature of the 
non-Gaussian covariance. The first (largest A) eigen- 
vector is the combination of bandpowers whose vari- 
ance degrades the most (relative to Gaussian) when non- 
Gaussianity is included; the second eigenvector degrades 
the second most, and so on. For {TT, TE, EE}, aU eigen- 
values are close to 1, and the variance degradation is 
essentially negligible, as quantified by the inequality in 
Eq. H28() . in any direction in parameter space. 

For BB bandpowers, we show the first few eigenval- 
ues and eigenvectors in Figure ^ The main effect of the 
non-Gaussian covariance is to degrade the variance, by a 
factor of ~ 10, for one KL component which is coherent 
across a wide range of I and has roughly the same shape 
as the fiducial BB spectrum. This is consistent with 
0, in which we found that non-Gaussianity degraded 
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FIG. 1: First three KL eigenmodes, defined by Eq. I27i . for the BB 
power spectrum. These represent modes in BB whose true variance 
is larger than the variance estimated from Gaussian statistics; the 
eigenvalue A is the ratio of the two. 



0.006 - 



0.004 



0.002 




500 1000 
Multipole 1 



1500 



FIG. 2: Principal components K\(r), K2{1) of the lensing poten- 
tial Cf^ obtained from CMB measurements to imax = 2000, as 
described in illV Al These represent modes in the lensing potential 
which are constrained by measuring cither lensed {T,E} or lensed 
B-modes respectively. 



the uncertainty in the overall amplitude of the BB spec- 
trum by a factor of ^ 10 when sample variance limited 
to Zniax = 2000. Here we see that this single statement 
roughly characterizes the entire non-Gaussian covariance 
between {TT, TE, EE, BB} bandpowers. 



IV. CMB LENSING OBSERVABLES 

In this section, we quantify the information in terms 
of the power spectrum of the lenses. An examination 
of the information contained in the power spectrum of 
the lenses serves a dual purpose: it is a parameterization 
independent quantification of the additional information 
from lensing ( ^lY Al IIV B|l . and it exposes the origin of 
the non-Gaussian covariance of the CMB as arising from 
the sample variance of the lenses (i JIV ( 
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A. Principal Components 

We begin by choosing the parameters of interest to be 
fluctuations pi in the power spectrum of the lenses around 
the fiducial model 

Cf^ = (f +pOCf*|fid. (34) 

With these parameters in the Fisher matrix of Eq. (|25|1 . 
the covariancc matrix 

(C)"' = (F-i)"' (35) 

can be interpreted as that of the measurements of Cf'^ 
under the assumption that the parameters that control 
the unlensed CMB are fixed. 

The principal components or eigenvectors of this co- 
variance matrix determine the best constrained linear 
combinations of Cf"^. We find that if lensing B- modes 
are not observed, the covariance C" is dominated by one 
well-constrained component Ki{l), which we show in Fig- 
ure |21 Equivalently, this means that only one observable 
in Cf"^ is constrained by lensed {T, E} power spectra: 

©^'"E 7^^1(0- (36) 

The power spectrum C^'^'^jfid of the fiducial model is 
scaled out of the weights Ki{l), such that deviations from 
01 = 1 represent the fractional change in the weighted 
amplitude of the power. Hence the normalization is cho- 
sen such that J2i -^i(0 = 1- 

On the other hand, if we make the artificial assumption 
that lensing i3-modes are observed but lensed {T, E} are 
not, then we find that the covariance is dominated by a 
single broad component K2, which peaks at / ~ 200 and 
includes a wide range of I. Lensed B-mode measurements 
therefore constrain a second observable 

Q^-Et^^^W- (37) 

The principal components Ki{l) were computed as- 
suming cosmic variance limited CMB measurements to 
^max = 2000; however, the shape of the eigenmodes re- 
mains nearly the same if Imax is lowered, or if a white 
noise power spectrum is used in place of a cutoff in 
I. Therefore, the observables Qi provide a parameter- 
independent representation of the information in the 
lensed CMB regardless of the noise characteristics. A 
caveat to this statement is that we never consider CMB 
multipoles beyond /max — 2000 in this paper; relaxing 
this assumption may permit additional modes in the lens- 
ing potential to be constrained. 

For BB, the higher principal components are not com- 
pletely negligible; the second-best constrained compo- 
nent has a variance which is worse than K2(l) by a fac- 
tor of 7. We have found that constraints from higher 




2 4 6 8 10 



z 

FIG. 3: Redshift sensitivity of the lensing observables ©i near the 
fiducial model. To good approximation the observables constrain 
the amplitude of cf "^ around multipoles near the median of the 
eignmodes of Fig. |^ Iki = 114, Ik2 = 440. The redshift sensitivi- 
ties Zi at these multipoles (see Eq. 1391 ) are plotted for the fiducial 
model. 

components can almost always be neglected in parame- 
ter forecasts, but can have some impact on degenerate 
directions involving curvature for a measurement of lens- 
ing i3-modes which is close to all-sky cosmic variance 
limited. In the rest of the paper, we will ignore higher 
components from lensed i3-modes. 

The structure of the two eigenmodes is related to the 
lensing kernels of Eq. Given that power in the de- 
flection angles peaks at h < 100 in the fiducial model, 
lensing mainly acts as a convolution kernel of width li 
on the high I CMB power spectrum. The {T, E} ker- 
nels share a similar structure since the angle between the 
lensed and unlensed 1 is of order li/l. The B kernel is 
weighted toward higher li for the same reason. Like- 
wise the dominance of a single mode in {T, E} refiects 
the tight range in li of the convolution compared with 
typical structure in the unlensed power spectra. 

B. Parameter Sensitivity 

Next to understand how sensitivity to these eigen- 
modes translates into cosmological parameters, let us ex- 
amine their construction in both the multipole and red- 
shift direction. The change in the observables due to cos- 
mological parameters can be derived from Eqs. H36|l - H37() 
once the change in Cl"^ is known. 

In Fig. 21 we plot the derivatives dCf"^ /dpa for several 
cosmological parameters Pa- The corresponding deriva- 
tives of the observables are given in Tab. ^1 Since the 
acoustic peaks constrain 

lA='7rDAiD,)/s,, (38) 

where Ss is the sound horizon at recombination Dg, we 
take these derivatives at fixed Ia (by adjusting J^de)- 
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They then quantify the additional sensitivity to cosmo- 
logical parameters introduced by lensing. 

Notice that the derivatives of the power spectra are 
quite flat compared across the multipoles where the two 
principal components have support (see Fig. Hence 
the sensitivity of the observables to most parameters 
can be accurately determined from the sensitivity of the 
power spectra at the median multipoles iKi of the princi- 
pal components, defined by X^'^i — 1/2: Iki — 114 
and Ik2 = 440. 

Next, to understand the relative sensitivities to differ- 
ent parameters, consider the fact that Cf"^ is determined 
by a projection of the matter power spectrum with a 
well-defined redshift sensitivity. In Fig. O we plot this 
sensitivity Zi{z), where 





Oi 


92 




-0.24 


-0.34 


d&i/dwo 


-0.14 


-0.12 


dQi/dWa 


-0.072 


-0.061 


dQ^/dQ.K 


-8.24 


-9.17 


de,/d{nch'^) 


17.0 


24.7 


96, /a In Sc 


2.00 


2.09 


de,/dlnlA 


2.37 


2.99 



TABLE II: Derivatives of the observables ©i, 02 with respect to 
parameters of interest (top) and nuisance parameters (bottom). In 
all rows except the last, the derivatives are taken adjusting Que to 
hold I A fixed. Units for Yl tt^u are eV. 



cf 



Cf 



fid 



dzZi{z) . 



(39) 



These weights are calculated under the Limber approxi- 
mation (see e.g. 0|). In the fiducial model / dzZi = 1, 
so that fluctuations in Zi determine fluctuations in the 
observables as 



^e,: « — = / dz6Z,{z). 



I fid 



(40) 



We expect this to be a reasonable approximation since 
5Cf'^ is not a rapidly varying function of / (see Figure0J. 

To make the above considerations more concrete, con- 
sider the sensitivity to changes in the distance Da, ex- 
pansion rate iJ, growth rate of the gravitational po- 
tential G and the shape of the matter power spectrum 
k^Pik) /27r3 at the lens redshift 



5Z^ 
Z, 



SD 



drf 3_ 



dlnAi 



SG 



dDA{D,-D) 



Da{Ds-D) 



d\nk 



k={lK./DA) 



(41) 



A typical value for the slope of the power spectrum gives 
rii ~ 1. 

CMB lensing is sensitive mainly to high redshift 
changes in the amount of lensing and correspondingly 
01 has a median redshift of z ~ 2 and Q2, z ~ 3. The 
observables best constrain cosmological parameters that 
change the geometry and growth in a coherent fashion at 
these redshifts. For example, a negative spatial curvature 
decreases the angular diameter distance to the lens red- 
shifts given the fixed and also decreases the growth 
rate. Both of these effects persist to 2; ~ 2 — 3 and hence 
lensing is highly sensitive to the curvature. Likewise mas- 
sive neutrinos slow the growth from the time they become 
non-relativistic near recombination. They also create a 
substantial change at high redshift. 

For the dark energy, we take its equation of state to be 
parameterized by {wq, Wa} according to Eq. H24|) . These 
parameters suffer in sensitivity in that the changes they 
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with respect to the parameters ^ nii, , 
wo, Wa, and Qk, illustrating the different I dependence. As in 
Tab, im the derivatives are taken adjusting Que to hold Ij^ fixed. 



FIG. 4: Derivatives of C, 



induce on the geometric and growth parameters are dom- 
inant at z < 1. Furthermore their effects on the expan- 
sion rate, distance and growth tend to cancel in the ob- 
servables (see [ijjl for a more extended discussion). Note 
that this a feature of the specific parameter set chosen 
and may not apply to all dark energy models. The princi- 
pal component technique allows a parameter independent 
way of quantifying the lensing information. 

The geometric and growth parameters are not the only 
ones that affect the observables. Since the observables 
are directly related to the power spectrum of the lenses 
any parameter that alters it also alters the observables 
causing parameter degeneracies (see ^V\ . The most im- 
portant of these are flch^ and the amplitude of the initial 
power spectrum (5^ . The sensitivity of the observables to 
these parameters are given in Tab. |nj For completeness, 
we also give the sensitivity to I a, though errors on this 
parameter from the acoustic peaks will be negligible for 
lensing purposes. 
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FIG. 5: Uncertainty on 0i from lensed T alone (top) and lonsed 
{T,E} (middle), and uncertainty on ©2 from lensed B (bottom), 
for varying beam size and noise level Ap. We assume Ay is given 
by Ap/\/2 throughout, including the top panel. Only multipoles 
below luiax = 2000 are included. For the zero beam cases, we also 
show the uncertainties that would be obtained if Gaussian statistics 
were falsely assumed (dashed). The impact of non-Gaussian contri- 
butions is negligible in {T, E} but significant for B. The horizontal 
lines are sample variance limits given by Eq. 1421 . 



C. Sample and Noise Variance 



The sample and noise variance of the lensing observ- 
ables determines the errors on cosmological parameters. 
Furthermore, this characterization of the errors illumi- 
nates the origin of the non-Gaussian bandpower covari- 
ance. 

In Figure |31 we show the Icr errors on Qi which are ob- 
tained for different instrumental sensitivities. Combined 
with the sensitivity of Qi to cosmological parameters, 
these results may be used to forecast parameter errors 
(see 

As the instrumental noise goes to zero, the errors satu- 
rate at the combined sample variance limits of the lenses 
and unlensed CMB. To understand the relationship be- 
tween sample variance and non-Gaussian covariance, con- 
sider first an ideal noise-free direct measurement of the 
lensing potential (p. From the definition (|36|l . the sam- 
ple variance limit on a measurement of the observables 



arising from the lenses is given by 

^iv(e.) = ^if.(Oi^.(/')(Csv)"' 



1,1' 



(42) 



Here we have taken Cgy 



2/skv'5ii'/(2^ + 1) under the 
assumption that (/) is a Gaussian field. 

In Figure 13 the sample variance limits of Eq. H42|) are 
shown as dashed horizontal lines. In the limit of zero 
noise, lensed CMB measurements can measure Qi and 
&2 at nearly their sample variance limits. Direct mea- 
surements of the <j) field, e.g. from CMB lens reconstruc- 
tion, do however contain more information than the ob- 
servables Qi . They are only two of many quantities that 
can be constructed with a sample variance limited Cgy. 
These other quantities are related to changes in the shape 
of Cf"^ and may be useful for breaking degeneracies. 

The two observables do however contain the majority 
of the low I information on the amplitude of Cf"^. To 



see this, consider the mode pi - 
parameter can be measured to 



A - 1 in Eq. This 



E 



21 + 1 



/sky 



(43) 



The qualitative difference between Eqs. H42II and H43|) is 
that the former is limited by the multipoles with the 
largest weighted sample variance. The latter is limited 
by the smallest and hence determined by the cutoff l^max- 
In the fiducial model (Tgy(A) = crgv(0i) at Z^max = 198 
and cr|v(^) = crivCQz) at Z^max = 705. Given that CMB 
power spectrum measurements out to /max = 2000 cap- 
ture nearly all of the information on the two observables, 
they also capture essentially all of the information on the 
amplitude of the lens power spectra near the median of 
the weights of the eigenmodes. 

Finally, if Gaussian statistics are falsely assumed for 
BB, then one would conclude that a high sensitivity mea- 
surement of CMB B-modes constrains the lensing observ- 
able Q2 better than the sample variance limit (see Fig.[SJ. 
From this perspective, one can gain intuition into why 
non-Gaussianity is significant for lensed B-modes. If the 
B-modes were perfectly Gaussian, then the overall am- 
plitude of the lensing could be constrained to within the 
sample variance of the smallest scale fluctuations in the 
CMB. In reality, the amplitude of the lensing is limited 
by the sample variance of the lenses near the median of 
the eigenvectors, i.e. on degree scales. Band powers be- 
low this scale covary in amplitude since the induced B- 
modes share the same lens fluctuations. This covariance 
becomes noticeable when the intrinsic sample variance of 
the ii^-modes becomes subdominant due to binning. 

The principal component analysis also shows why non- 
Gaussianity is not a significant limitation for {T, E}. The 
Gaussian errors on Oi from observations of lensed {T, E} 



9 





V 


Sfwhm 


At 


Ap 


/sky 


Planck 


100 GHz 


9.2' 


51 




0.8 




143 GHz 


7.1' 


43 


78 


0.8 




217 GHz 


5.0' 


65 


135 


0.8 


Deep 10% 




1.0' 


1.00 


1.41 


0.1 



TABLE III: Assumed experimental specifications for the reference 
survey in Jvl consisting of Planck and a deep polarization field with 
10% of sky. The noise parameters Ay and Ap are given in units 
of /iK-arcmin. For the combined reference survey, (t(©i) = 0.025 
and o-(e2) = 0.008. 



(Fig. 13 never exceed the sample variance limit. This is 
because the sample variance of the high I unlensed {T, E] 
fields still dominate the measurement of the small frac- 
tional changes induced by lensing. The B field does not 
suffer from this problem in that at high I it is completely 
generated by lensing. 



V. PARAMETER CASE STUDIES 

We now study parameter constraints from the lensed 
CMB, with non-Gaussian contributions to the power 
spectrum covariance included. In general, parameters 
which affect the high-redshift universe or the angular di- 
ameter distance to last scattering are well measured even 
without lensing. Lensing mainly helps in breaking degen- 
eracies that leave the observables at recombination fixed. 
We illustrate this degeneracy breaking with massive neu- 
trinos Tny\^ a constant dark energy equation of state 
{w}, equation of state evolution {wojWa}, and spatial 
curvature \^k\. Given that the lensed CMB adds two 
new observables as described in the previous section, we 
will study these additional parameters two at a time. We 
start by examining the rtiy^ w\ case extensively, as an 
illustration of the impact of non-Gaussianity, as well as a 
worked example of the use of the lensing observables 8^. 

Throughout this section, we will consider a reference 
survey fTab. llTT|) which consists of Planck measurements 
on most of the sky, combined with a smaller deep survey 
which measures lensing of the temperature and polariza- 
tion with higher signal-to-noise. (We assume that the 
deep survey area is a subset of the Planck survey area.) 



A. Neutrinos and constant equation of state 

We consider here the constraints in the {^r7ii,,w} 
plane at fixed employing a direct Fisher matrix fore- 
cast in the parameter space and through the intermediary 
measurement of the lensing observables, first by adding 
its information to the unlensed spectra in the full pa- 
rameter space and then by pre-marginalizing parameters 
that control the high redshift physics at recombination. 



1. Direct forecasts 

In Figure |B| left panels, we show the constraints 
falsely assuming Gaussian statistics. For illustrative 
purposes, we first show the errors with the parameters 
{i^ch^, In S(^}, fixed to their fiducial values (top left), but 
with the remaining parameters {r^oE, f^fe/i^, t, ris, r} 
marginalized. As we shall see, the former two parame- 
ters affect the lensing observables and cause degeneracies 
with the parameters of interest. 

The large surrounding ellipse in the top left panel rep- 
resents the double angular diameter distance degener- 
acy expected from the unlensed CMB. Adding informa- 
tion from lensed {T,E} (horizontally shaded) constrains 
a slightly different combination of degenerate parame- 
ters than adding information from lensed B (vertically 
shaded). This is because changes in the two parameters 
{^m^,-u;} affect Cf"^ differently in I (Fig.0J and hence 
the two lensing observables described in HIVI 

We next consider the effect of including the non- 
Gaussian covariance, in the top right panel of Fig. 
Compared to Gaussian uncertainties, the constraints 
from lensed {T, E} are unaffected (as expected from 
^III B|) . while constraints from lensed B are degraded 
for the best constrained combination of the two param- 
eters. Note that because we have considered a nearly 
degenerate two parameter space, even here the effect of 
non-Gaussianity is hidden from the errors on a single pa- 
rameter marginalized over the other. 

However, the discussion so far has assumed perfect 
priors for the parameters {flch^ 6(^}. If these are 
marginalized with only information from the reference 
survey, the effect of non-Gaussianity is overwhelmed by 
the effect of marginalizing (see Fig. |H| bottom panels). 
This is true even if B were the only source of lensing in- 
formation or if there were only one additional parameter. 

Nonetheless the effect of non-Gaussianity in the full 
parameter space is not negligible; it does enlarge the vol- 
ume which is allowed. The enlargement occurs along a 
direction which is a combination of several parameters 
{r2c^^,ln (5^, w, ^ TOiy}. Though hidden by marginaliza- 
tion, this degradation can be exposed if there are exter- 
nal prior measurements of other combinations of these 
parameters. It can also reduce the apparent goodness of 
fit to the best-fit model (see Appendix . Finally for 
different reference surveys, e.g. a smaller but deeper sur- 
vey, a Gaussian assumption can again give a misleading 
answer on the utility of the _B-modes (see Will . 

These results can be reproduced and understood using 
the lensing observables Qi as we shall now see. 



2. Additive observables forecasts 

The lensing observables provide a general framework 
for forecasting the additional constraints supplied by 
lensing in any parameter space. We can use the direct 
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FIG. 6: Degeneracy breaking in the —w plane from CMB lensing, for the reference survey in Tab. lIIll Ellipses here and throughout 

are plotted at Ax'^ = 1 and not 68% CL. In each panel, the surrounding ellipse represents parameter constraints from unlensed {T, B}, 
the blue/horizontally-shaded ellipse represents constraints from lensed {T, E}, the red/ vertically-shaded ellipse represents constraints from 
unlensed {T, E} + lensed B, and the inner solid ellipse represents constraints from lensed {T, E, _B}. In the top panels, we show Gaussian 
(top left) and non-Gaussian (top right) constraints with {f!c/i^, In 6^} held fixed, and the remaining parameters {Q,f,h? ,Q.Y)^,T,ns,r} 
marginalized. In the bottom panels, we show Gaussian (bottom left) and non-Gaussian (bottom right) constraints with all parameters 
marginalized. 




-0.4 -0.2 0.2 0.4 0.6 
Em, (eV) 

FIG. 7: Forecasted errors in {^m^,w} as in Fig. l6l flower right 
panel, non-Gaussian, fully marginalized) but with the additive lens- 
ing observables prescription from iiV A 21 Note the excellent agree- 
ment between the direct and observables approaches. 



forecasts from the preceding section as a basis for com- 
parison with the foUowing construction: 

1. From the parameters of the survey under consid- 
eration (sky coverage, noise, and beam), compute 
uncertainties (see Fig.O on the observables Qi and 
assume they are independent, obtaining a 2x2 co- 
variance matrix: 

Co„e.,e,). J. («, 

2. Transform these parameter errors into a Fisher ma- 
trix in the desired parameter space with 

F':^' = (5„e')Cov(e„ Qj)-\dpQ^) (45) 

and add this to the Fisher matrix of the unlensed 
CMB or any external data set. 
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Some of the more common parameter derivatives for 
use in Eq. H45|l are given in Tab. ^ Others can be ap- 
proximated by evaluating Cf'^ at the median multipole. 
An advantage of the observables scheme is that it auto- 
matically includes the effect of non-Gaussianity without 
reference to the bandpower covariance matrix. One sim- 
ply uses the non-Gaussian errors from Fig. 

The main caveat to the observables prescription is that 
it implicitly assumes that lensing is an additive source of 
information. It should not be used to forecast parameters 
for which lensing destroys information such as tensors or 
features in the initial power spectrum. In the particular 
case of the tensor-to-scalar ratio r, we find that non- 
Gaussianity in the lensed CMB is always negligible when 
forecasting uncertainties, even when low-redshift param- 
eters are marginalized. 

We show the w — example in Fig. [7| Note the 

excellent agreement with the direct Fisher calculation in 
Fig. (lower right panel) in all respects. 



3. Nuisance-marginalized observables 

This prescription can be further simplified so as not to 
require computation of the unlensed CMB Fisher matrix 
or external information for forecasting parameter errors 
which rely mainly on the lensing information. In this case 
the extra information is utilized to remove the parameter 
degeneracy from f2c^^ and In S(; in the observables O^. 
These parameters can be premarginalized and dropped 
from the Fisher matrix constructed from O^. 

Operationally, to pre-marginalize the nuisance param- 
eters, we add the "nuisance errors" 



C0Vnuis(6i,ej 



^ ^^^-.^(in.c) 



9 In 9 In 5, 



C 



to the "measurement errors" defined by Eq. (|44|l . to get 
an effective covariance matrix Covofr — (Gov + Govnuis). 
We then use Gov^ff instead of Gov in Eq. H45|) . to obtain 
a Fisher matrix in a parameter space where {ilch^, In 5^} 
is excluded. 

In (|46|l . we have neglected the correlations between 
ilch^ and In Sq but these provide a negligible effect if 
the source of information is internal to the GMB. Note 
that Govcff will typically show a high degree of corre- 
lation between the two observables since the increase in 
the effective errors is along directions that are degenerate 
with the two parameters. For parameters that do not in- 
duce a degenerate change, this correlation reflects extra 
information (see WI B|) . 

Let us illustrate this pre-marginalization scheme for 
{'^m^,w} and the reference survey. For comparison 
with Fig. El we will compute errors with priors on 
{rich^^hi (5^} that fix them completely, and priors that 
are derived from determinations internal to the reference 
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0.01 
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0.18 
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O"nuis(02) 


0.25 
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Corrnuis(6i, 02) 


0.99 


0.99 


0.99 


0.99 


CTcff (6l) 




0.054 


0.030 


0.013 


0-cfl(62) 






0.025 


0.013 


Correff(ei,e2) 






0.52 


0.71 



TABLE IV: Uncertainties in lensing observables from third-year 
WMAP Ifij, Planck alone fTab. lTTTl top), the reference experiment 
consisting of both parts of Tab. 11111 and an ideal survey which is 
all-sky cosmic variance limited in temperature and polarization to 
'max = 2000. Top: Measurement errors on the lensing observables 
{©1,02} computed from the raw sensitivity to lensed {T,E} and 
B respectively, as in Fig. |5] Middle: "Nuisance" errors of Eq. I46i 
on lensing observables, computed by propagating each survey's un- 
lensed uncertainties on {^^h^, In S^} (also shown) using derivatives 
from Tab. [H] Bottom: Total errors on lensing observables, com- 
puted by adding 2x2 covariance matrices {CoVgff = Cov-|-CoVnuis). 
As described in the text, these are the effective errors on lensing ob- 
servables when constraining low-redshift parameters with nuisance 
parameters {Oc?t^,ln <5^} premarginalized. 



survey (see Tab. IIVII . The former case also corresponds 
to using Gov (Eq. (03) in place of Govoff. Gombined 
with the parameter derivatives from Tab. ^ the result 
of this procedure is shown in Fig. |S1 

The constraints agree well with the direct Fisher cal- 
culation in the right panels of Fig. along the best con- 
strained direction. For example, with ^ fixed, the 
uncertainty a{w) is 0.150 from the direct Fisher calcu- 
lation ( ^V A ip . 0.170 from the additive-observables pre- 
scription f W A2p . and 0.178 from the premarginalized- 
observables prescription described in this section. In 
practice, the main difference between the last two is 
that in the premarginalized-observables prescription, the 
constraints computed are from lensing alone; informa- 
tion from the unlensed GMB which helps constrain 
intermediate-redshift parameters (such as ISW) is not 
included. 



A pedagogical advantage of pre-marginalization is that 
it clarifies the limiting source of uncertainty. For ex- 
ample, in our reference survey, o'nuis(02) is substantially 
larger than cr (82). The utility of the lensing observable 
is therefore limited by the nuisance parameters and not 
by the sample or noise variance on lensing. This also 
explains the negligible impact of non-Gaussian sample 
variance on marginalized parameter errors in Fig. 
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FIG. 8; Forecasted errors in m^, ui} as in Fig. lol fright panels, 
non-Gaussian) but with the premarginalized-observables prescrip- 
tion from i|V A 31 The blue/horizontally-shadod regions show only 
the 01 constraint from lensed {T, i?}, the red/vertically-shaded 
regions show only the O2 constraint from lensed B, and the solid 
regions show the combined constraint from both observables. The 
top panel shows the case where {Q,c.h? ,\n (5^} are completely fixed 
by the external prior; the bottom panel shows the case where the 
parameters are internally determined by the reference survey itself. 



B. Evolution of equation of state 

We now consider CMB constraints on a time- 
dependent dark energy equation of state through 
{wo,Wa} [see Eq. (1^ ]. For purposes of this subsection, 
we assume that both the neutrino mass and spatial cur- 
vature are fixed and marginahze {Ode, i^bh"^, ^ch^, t, 
ris, r}. 

In Fig. 13 we show errors on wq and Wa, for the ref- 
erence survey of Tab. IIIII using the premarginalized ob- 
servables scheme from W A3I In the top panel we show 
the direct Fisher matrix calculation for comparison. One 
parameter combination is constrained by CMB lensing, 
but the complementary direction is degenerate, acquiring 
CMB constraints only from the unlensed CMB through 
the ISW effect (c.f. |l2). 

The {wojWa} degeneracy can be seen directly by not- 
ing that the wq and Wa derivatives, taken at constant 
angular diameter distance, of Cf"^ are nearly propor- 
tional (Fig. 01) , so that a parameter space direction exists 
which preserves both Cf"^ and the unlensed CMB. This 
is also seen in Tab. m where the wq and Wa derivatives 
of {6i, 62} are nearly proportional. 




FIG. 9: Joint uncertainties on dark energy parameters {wa,wo}, 
for the reference survey in Tab. |m] with {Y^m^,^!^} fixed to 
their fiducial values and all other parameters marginalized. The 
top panel shows the direct Fisher matrix calculation; the bot- 
tom panel shows errors calculated using the premarginalized- 
observables scheme f il V A 31 . In this parameter space, the lensing 
observables are degenerate; the ©i constraints (blue/horizontally- 
shaded), the 02 constraints (red/vertically-shaded), and the com- 
bined constraints (solid) are nearly identical. Given this nearly 
perfect degeneracy, the Fisher forecast should not be interpreted 
literally along the full extent of the degeneracy. 



The degeneracy can be interpreted as the statement 
that the lensed CMB constrains w{z) mainly around the 
pivot redshift Zp ~ 1 with uncertainty a{wp) = 0.15 for 
the reference survey. The pivot redshift can be inter- 
preted as a combination of the high redshift weight of 
CMB lensing discussed in ^IIVBI and the dark energy 
parameters which change observables strongly only at 
z < 1. Note that this pivot redshift is substantially 
higher than many other cosmological probes and im- 
plies that CMB lensing will provide complementary con- 
straints on the evolution of the equation of state. We 
also note that the pivot redshift and a{wp) will depend 
on the underlying fiducial model. 



C. Curvature 

Given the sensitivity of CMB lensing to changes at 
high redshift, constraints on spatial curvature are much 
stronger than those on the dark energy when compared 
with low-redshift probes of the expansion history. In 
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Fig- El we show the reference constraints in the {^k, w} 
plane with {fli,h^ ,rtch^ ,rti:i^,T,ns,r} marginahzed and 
{J2 m^yWa} fixed. 

The two observables are again nearly degenerate in this 
plane. The direction of the degeneracy is well but not 
perfectly matched in the observables scheme for the case 
of 81. The reason for this is that the unlensed CMB 
carries information on the curvature both from the ISW 
effect and from the intrinsic sharpness of the acoustic 
peaks (see [23; Fig. 11). The latter effect comes from 
the geometrical projection of fc-space power to Z-space 
through the radial eigenfunctions. Unlike lensing this ef- 
fect smooths (negative curvature) or sharpens (positive 
curvature) even the low order peaks as it is associated 
with the curvature across the last scattering surface as a 
whole. For the high order peaks it destructively interferes 
with the lensing effect and mildly violates the assump- 
tion that lensing is an independent additive source of 
information. The interference and the rotation it causes 
however are small and, if desired, can be accounted for 
by a ^ 20% lowering of dQi/dil.K for surveys that utilize 
information near ^max = 2000. 

More importantly, the direction that is well con- 
strained has a large fix and only a small w component, 
i.e. the degeneracy is very steep in the curvature direc- 
tion. The implication is that it takes only a weak exter- 
nal constraint on w to break this degeneracy completely. 
When combined with other dark energy probes, the lens- 
ing observables can be thought of as fixing the curvature. 
We explore this use of CMB lensing further in a separate 
piece |2ll| . 

Even given other probes that break the {^Ik,w} de- 
generacy, fix remains nearly degenerate with ^TOi^. In 
Fig. ^2 we show constraints in this plane with {wq, Wa} 
fixed. The observables approach again accurately mod- 
els the well constrained direction aside from the slight 
rotation of the @i constraint which makes curvature 
nearly perfectly degenerate with neutrinos. Breaking this 
degeneracy externally will require independent probes 
with limits of (j(CIk) ^ 0.005 (perhaps with measure- 
ments of the angular diameter distance to 2: ~ 3 [2^ ) or 
o-(E™^) < 0-2eV. 



VI. SURVEY OPTIMIZATION 

In the previous section, we quantified the information 
supplied by CMB lensing in the context of a specific ref- 
erence survey. We conclude this paper by considering 
how surveys can optimize the extraction of cosmological 
information from lensed power spectra. 



A. Qch'^, r, and /max 

We have seen in ijlVI that imperfect knowledge of 
{flcfi^, In S(^} results in nuisance errors which limit the 




FIG. 10: Joint uncertainties on {Uk,wo}, for the reference sur- 
vey of Tab. mil with {Y^ m,^,Wa} fixed to their fiducial values and 
all other parameters marginalized. The top panel shows the direct 
Fisher matrix calculation; the bottom panel shows errors calculated 
using the premarginalized-observables scheme from i|V A 31 As ex- 
pected, the two agree in the well constrained direction, but the 
second scheme overestimates the errors in the poorly constrained, 
ISW-limited, direction. 



ability to extract cosmological information from lens- 
ing. The nuisance errors anms{Oi) represent a floor be- 
low which improving sensitivity to each lensing observ- 
able individually does not improve cosmological param- 
eter constraints. More precisely, improving <7{Qi) be- 
yond this level will correlate low-redshift parameters to 
{flch^, In (5<^}, in such a way that marginalized uncertain- 
ties do not improve if only one of the observables is mea- 
sured. Joint measurement can slightly improve on this 
"floor" due to the difference in how the nuisance param- 
eters affect the two. 

One source of nuisance error arises from uncertainty 
in flch^ as measured from the primary CMB or exter- 
nal data. This is the dominant source of error through- 
out Tab. IIVI in which we have assumed that the range 
of multipolcs from which cosmology can be extracted is 
limited to /max = 2000. This cut reflects an estimate of 
possible contamination from other secondaries and fore- 
grounds and is currently uncertain. Fortunately the de- 
pendence of a{flch^) on /max is not particularly sensitive 
strong (see Fig. IT^ , provided that it exceeds the knee 
at /max ~ 700, corresponding to the trough between the 
second and third acoustic peaks. 

Another conclusion from Fig. E|is that achieving cos- 
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FIG. 11; Joint uncertainties on dark energy parameters 
i^i^^^k} with {wo,Wa} fixed to their fiducial values, for the 
reference survey of Tab. 11111 using the direct Fisher calculation 
(top panel) and the premarginalized-observables scheme (bottom 
panel). (In the top panel, the red/ vertically-shaded ellipse repre- 
senting unlensed {T, E} and lensed B is hidden by the solid ellipse 
representing lensed {T, E, B}.) 
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FIG. 12: Uncertainty <j{flch'^) from CMB measurements with 
varying imax, for Planck sensitivity (Tab. ITTTl . and ideal measure- 
ment of unlensed CMB temperature and polarization. 



mic variance limited measurement of i?-mode acoustic 
peaks at intermediate I would significantly improve flch^ 
constraints from Planck. For example, an ideal mea- 
surement of {T, E} with Zmax ^ 500 would obtain flch^ 
constraints comparable to Planck, even though Planck 
will measure CMB temperature anisotropy well into the 



damping tail. We have found that the inability to extract 
flch^ constraints from the temperature damping tail is 
due to confusion with the spectral index Us ; if Ug is fixed 
rather than marginalized, then Planck's flch^ constraints 
significantly improve. 

As Zinax is increased in Fig. ^1 the error on flch^ im- 
proves, but the error on our second nuisance parameter 
In stays nearly constant. This is because the degener- 
acy between In and r in the unlensed CMB is broken 
only by the reionization signal at low polarization multi- 
poles. Even at Zmax = 2000, the largest value we consider 
in this paper, the nuisance errors are dominated by un- 
certainty in Q,ch?. This may suggest that In 5q is never 
important as a nuisance parameter. There are three rea- 
sons why this may not be the case in practice. 

First, errors on Vlch^ can be in principle be improved 
over CMB determinations by external sources such as 
weak lensing of galaxies or even from the CMB itself 
through lens reconstruction, whereas amplitude uncer- 
tainties from reionization will likely remain. Second, we 
assume that large angle polarization foregrounds can be 
perfectly removed to the cosmic variance limit through 
the frequency channels not used for cosmology. These 
foregrounds are now known to dominate the polarization 
signal in all bands 0|. 

Third, our parameter forecasts so far have implicitly 
assumed sharp reionization, characterized by the sin- 
gle parameter r, the total optical depth to recombina- 
tion. More general models of reionization can include 
additional parameters which degrade uncertainties on r 
[2^ 01 J and therefore our nuisance parameter In 5q, be- 
yond what is shown in Tab. IIVI 

If we conservatively assume that the ionization his- 
tory takes an arbitrary form, uncertainties rise to <j{t) = 
cr(ln 5(^) = 0.01 and by Tab. m the corresponding 
nuisance errors are crnuis(0i) = o'nuis(02) = 0.02. Com- 
paring with Tab. IIVI this would be a comparable source 
of nuisance error to Vlch? for Planck, and would represent 
the dominant uncertainty for an ideal experiment that 
achieves cosmic variance on the i?-modes (c.f. Fig. 1121) . 

Finally the choice of /max also affects constraints on 
01. For example scaling back to /max = 1500 would 
degrade (t(8i) by 1.2 for Planck and 1.6 for the ideal 
measurement. 



B. Optimizing sensitivity to lensing B- modes 

An important issue for upcoming polarization exper- 
iments is optimizing sky coverage with the total inte- 
grated sensitivity fixed. In general this is not a well-posed 
question since the sky coverage can be optimized with 
respect to systematic errors, foreground contamination, 
sensitivity to the i?-mode power spectrum, sensitivity to 
tensor i3-modes, detection of lensing B-modes, i3-mode 
reconstruction of the lensing fields, or sensitivity to cos- 
mological parameters from the B-mode power spectrum. 
Here, we consider only the last of these, with particular 
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fixed integrated sensitivity ^ = 2 nK. The solid line shows the 
result using non-Gaussian statistics for lensed B-modes f jlTl : the 
dashed line shows the result if Gaussian statistics are falsely as- 
sumed. The horizontal lines represent the "nuisance" errors on ©2 
from imperfect measurement of {Qch^, In S(^}, for Planck and for 
ideal measurement of the unlensed CMB (Tab. ITVt . 



attention to how optimizing the sensitivity is affected by 
non-Gaussian statistics. 

The integrated sensitivity fi has units of temperature 
and is given by v/V NT where v is the instantaneous sen- 
sitivity (in mK-^sec) per Stokes parameter per detector, 

is the number of detectors, and T is the total integra- 
tion time. The noise variance per steradian in Eq. 115|) 
is then given by Ap = fi\/A where A is the survey area 
in steradians. 

In ijlVI we found that aU cosmological constraints from 
the lensing B-mode power spectrum are derived from the 
single observable 02. Therefore, there is a natural figure 
of merit for optimizing sky coverage: the la error (t(02). 
In Fig. 1131 we show (t(02) for varying /sky, assuming 
fixed integrated sensitivity /i = 2 nK. Incorporating non- 
Gaussian statistics increases the optimal value of /sky 
by a factor of 3 relative to Gaussian, and significantly 
steepens the dependence of (7(82) on /sky in a manner 
which disfavors small /sky 

We have shown the /sky optimization in detail for inte- 
grated sensitivity fi — 2 nK; in general, the optimal /sky 
will scale with /i as /sky oc The same scaling is ob- 
tained assuming either Gaussian or non-Gaussian statis- 
tics; therefore, the optimal patch size with non-Gaussian 
statistics incorporated is a factor of 3 larger than the 
Gaussian value, independent of the integrated sensitiv- 
ity. Another way of stating the optimality criterion is: 
the optimal patch size is always chosen so that the noise 
per steradian takes the value 



(Ap)optimai = 4.7 ^K-arcmin. 



(47) 



This criterion makes no reference to the value of /i but 
does assume zero beam; we have found that (Ap)optimai is 
nearly independent of beam size, provided that ^fwhm < 
15 arcmin. (We note that if Gaussian statistics were 
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Reference 
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a{wo) 


0.31 


0.18 


0.15 


0.07 




0.65 


0.38 


0.30 


0.15 


cr(nx) 


0.0076 


0.0032 


0.0025 


0.0013 




0.20 


0.085 


0.063 


0.032 



TABLE V: Uncertainties on each of {uiq, uij,, f^K, 52 sepa- 
rately with the others fixed, and high-redshift parameters marginal- 
ized. Here, "Deep5%" stands for a deep survey with /sky = 0.05, 
zero beam, and Ap = \/2At = 4.7 /iK-arcmin (see Eq. 1471 '). The 
reference survey is as in Eland covers /sky = 0.10, and "Ideal" 
refers to all-sky cosmic variance limited {T, E, B} to imax = 2000. 
Forecasts here are from the direct Fisher approach. 



falsely assumed for lensing _B-modes, then one would ob- 



tain (Ap 



optimal 



2.8 /xK- arcmin.) 



Another conclusion of the previous section was that, 
for sufficiently precise measurements of lensing i?-modes, 
the ability to extract cosmological information from the 
CMB alone is primarily limited by uncertainty in flch^ 
from the primary CMB and secondarily limited by reion- 
ization if the ionization history is complex. More pre- 
cisely, when the measurement error a{&2) becomes as 
good as the the nuisance error crnuis(02), then improved 
sensitivity to lensing B-modes serves mainly to correlate 
low-redshift and nuisance parameters, rather than im- 
proving marginalized uncertainties on either. In Fig. 1141 
we have shown the dependence of (t(02) on total sensi- 
tivity and sky coverage. 

We now consider optimization of a deep ground-based 
polarization survey designed to complement Planck. For 
the Planck prior on il.ch^, we found (Tnuis(02) = 0.028 
(Tab. IIVI) . Comparing with Fig. ^1 it is seen that a 
narrow-beam polarization survey with integrated sensi- 
tivity fi— 1-2 nK, and covering a few percent of the sky, 
will achieve (t(02) « crnuis(02) and is therefore nearly op- 
timal for extracting cosmological information from lens- 
ing i?-modes alone, within the limits of the Planck prior. 

This conclusion that a few percent of the sky is optimal 
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FIG. 15: Uncertainty on the single low-redshift parameter w, 
marginalized over high-redshift parameters, for Planck comple- 
mented by a deep survey with varying sky coverage /sky For the 
deep survey, we assume At = Ap/\/2 and consider three sensitiv- 
ity levels as indicated. 



for extracting the B-mode information alone is slightly 
modified once the T and E information from Planck and 
the deep ground-based survey itself are considered. The 
T and E measurements supply Qi and the nuisance pa- 
rameters. The first modification is that if Planck suc- 
ceeds in measuring Qi to the forecasted a{Qi) — 0.05 
then it alone has lensing information that is comparable 
to an /sky ^ 0.01 S-mode survey. The second mod- 
ification is that as the deep survey improves the nui- 
sance parameters directly or through 0i, the nuisance 
floor on the _B-modes correspondingly drops. Both con- 
siderations favor slightly larger /sky. To illustrate this 
we show in Fig. El the direct forecast on a single low- 
redshift parameter, e.g. w, if Planck is complemented by 
a deep survey with varying /sky For /sky less than the 
"knee" at /sky ^ 0.05, the uncertainty in w from the 

— 1/2 

-B-mode measurement itself improves as /sky but only 
becomes stronger than the lensing constraint from Planck 
for /sky ^ 0.01. For /sky ^ 0.05 improvements scale more 

— 1/3 

slowly as ~ /^j^^' . Half the total improvement comes 
from /sky < 0.1. In fact, a deep survey with /sky = 0.05 
can improve Planck lensing uncertainties on any one of 
{wq, Wa, ^ TOiy}, with the others fixed, by a factor 
of ~ 2 (Tab. 0) through measurement of the _B-mode 
observable. Moreover if lensing constraints from Planck 
prove impossible to extract due to foreground and sec- 
ondary contamination this improvement represents an- 
other factor of ~ 2 in errors. 



VII. DISCUSSION 

We have provided a comprehensive study of the addi- 
tional cosmological information supplied by lensed power 
spectra of the CMB temperature and polarization fields 
including the non-Gaussian covariance between band- 



power estimates. This covariance originates from the 
sample variance of the degree scale lenses on the CMB 
fields at smaller scales. It is nearly irrelevant for the tem- 
perature and ^^-polarization fields out to Z„iax = 2000 due 
to the larger sample variance of the unlensed CMB. For 
the amplitude of the _B-polarization field, it increases the 
variance by up to a factor of ~ 10 and changes the opti- 
mal observing strategy to one that covers a factor of ~ 3 
times more sky area. 

The impact of non-Gaussianity on parameter estima- 
tion as well as the net information content of the lensed 
spectra is more subtle. These answers depend on the 
choice of parameters and the external priors associated 
with them. We have provided a framework of lensing 
observables that greatly simplifies these examinations. 

In this framework, lensed CMB power spectra provide 
information on only two observables, one which deter- 
mines the lens power spectra at / ^ 100 associated with 
the {T, iJ}-fields and one which determines it at Z ^ 500 
associated with the _B-field. The observables are con- 
structed from the principal components of the lensing 
power spectrum Cf"^ . Non-Gaussianity is then automat- 
ically incorporated in the errors on the observables which 
will eventually approach, but never exceed, the sampling 
errors of the lenses as the measurements improve. 

This construction also illuminates the origin of param- 
eter degeneracies which can rapidly become the limiting 
source of uncertainties for parameters of interest. Any 
combination of parameters that leaves the lensing ob- 
servables and the CMB at recombination fixed within the 
errors cannot be determined. To illustrate these effects, 
we have isolated two parameters ilch^ and In that de- 
termine the shape and amplitude of the matter power 
spectrum respectively, and marginalized their uncertain- 
ties assuming internal CMB determinations of each from 
the Planck satellite. These become the limiting uncer- 
tainties once the observables are determined to the sev- 
eral percent level and are only slowly improved as the 
lensing survey itself improves the nuisance errors. While 
^Ich^ constraints can be improved externally to the CMB, 
those on In 5^ are more difficult to improve and may be 
limited by our understanding of reionization. 

There are also degeneracies within the space of the 
parameters of interest that control the expansion rate 
and growth of structure at intermediate redshifts. When 
taken one at a time, uncertainties on the parameters 
{'Wo,Wa,i^K,'^'mu} can be improved by a factor of ~ 
2—3, relative to Planck alone, by a deep ground-based po- 
larization survey on 5-10% of the sky. However {wq, Wa} 
are nearly perfectly degenerate in the lensing observables 
as are {r2i^,^m^} separately. The degeneracy between 
two parameters in each pair is weakly broken by the 
two observables. For example, when errors on ^ rrii, are 
marginalized over wq they degrade by a factor of 2 for the 
reference survey (see Fig. (H)). However sensitivity to the 
{^K, X] "^i'} PS-ir is much greater than to the dark energy 
parameters due to the high redshift weights of the lensing 
observables. When combining lensed CMB power spec- 
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tra with other more incisive probes of the dark energy, 
lensing essentially fixes one well-defined combination of 

Our conclusions have several caveats associated with 
them. The observables framework implicitly assumes 
that lensing is an independent and additive source of cos- 
mological information that may be combined with the 
intrinsic CMB anisotropy. An important exception to 
this statement occurs for tensor modes, where lensing 
B-modes mask the intrinsic _B-modes. Forecasts for ten- 
sor modes should be made employing lensed power spec- 
tra as a destructive contribution but here the Gaussian 
approximation suffices. The conversion between instru- 
mental noise and errors on the observables depends only 
mildly on the fiducial model given current cosmological 
constraints but we give a crude scaling in Appendix IbI 

Secondly, we have considered only the information con- 
tained in the lensed power spectrum. Beyond the power 
spectrum, non-Gaussianity from lensing allows a direct 
reconstruction of the lensing fields 0, |a| which carries 
substantially more information that can break parame- 
ter degeneracies 0, 0| . It may also allow "de- lensing" 
techniques that recover the intrinsic i?-modes from ten- 
sor modes [2^ . However techniques have yet to 
be developed that can remove systematics and contami- 
nation at the levels required. 

Thirdly, our parameter forecasts employ the Fisher 
matrix approximation. It is well known that Fisher ma- 
trix forecasts are not accurate along ill-constrained direc- 
tions in the parameter space. Hence our results are only 
robust for quantities that lensed power spectra constrain 
well. Finally, we never consider CMB multipoles beyond 
^max = 2000 in this paper. Well beyond this limit there is 
extra information on the high multipole structure of the 
lensing field but this is likely to prove difficult to extract 
in the presence of other secondaries and foregrounds. 
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APPENDIX A: GOODNESS OF FIT 

Given that lensing non-Gaussianity induces covariance 
in the band power estimates, it is interesting to ask 
whether the correct cosmological model would be in- 
ferred to be a bad fit to the observed bandpowers if 
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FIG. 16: Non-Gaussian corrections A(x^) (top panel) and 
AVar(x^) (bottom panel), as given by Eqs. |X3j, J^. We assume 
sample variance limited measurements of TT (black), EE (blue), 
and BB (red) up to maximum multipole imax- For fixed iniaxi the 
corrections are nearly independent of the number of bands; we il- 
lustrate this by showing a many-band fit (solid) and a one-band fit 
(dotted) for each power spectrum. The corrections are negligible 
for TT and EE but important for BB, if Ny^i^n^n is small and (max 
is large. 

non-Gaussian correlations were not included in the x^- 
This question was raised in the context of the first-year 
WMAP analysis Correspondingly, we define the 

naive statistic as a sum over bands, 

x' = E%/^> (Al) 

where the bandpower average (A^) is computed using 
lensed power spectra, and the Gaussian variance Qa is 
computed assuming lensed power spectra and Gaussian 
statistics, as in Eq. 

The statistic defined by Eq. (|A1|) fully incorporates 
the effects of lensing at the power spectrum level but 
neglects the non-Gaussian covariance between bandpow- 
ers. With non-Gaussianity included, the distribution is 
no longer a perfect x^, but acquires corrections 

(x') = ^dof + A(x'>, 
Var(x') = 2iVdof + AVar(x') , (A2) 

where the excess contributions A(x^), AVar(x^) arise 
only from higher-point correlations in the lensed CMB. 
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In this appendix, we study the size of these contributions, 
as a way to quantify the impact of non-Gaussianity. 

The first contribution in Eq. (jA2p can be written in 
terms of the bandpower covariance defined in Eq. (|12|) : 



(A3) 



In contrast, the full non-Gaussian contribution to 
Var(x^) is an eight-point correlation between CMB fields, 
and the results of this paper do not permit every term 
to be computed. However, if we make the approximation 
that the bandpowers are Gaussian variables, then it 
is given by: 



AVar(x^ 



2 ^ MjMj 



GiiG 



(A4) 



Since each bandpower is an average over many Fourier 
modes (see Eq. ©), the central limit theorem implies 
that this should be an accurate approximation. This gen- 
eral observation shows that, in the limit of wide bands, 
the bandpowers should always behave as Gaussian 
variables; lensing simply induces a Gaussian covariance 
between the bandpowers. 

We have found that the non-Gaussian contributions 
A(x^), AVar(x^) to unreduced values are nearly inde- 
pendent of the number of bands or degrees of freedom. 
In Figure [TEl we show these contributions for lensed TT, 
EE and BB power spectra, and for varying Zmax, in two 
extreme cases: a "many-band" fit with Al — 10, and 
a "one-band" fit across all multipoles up to /,nax- The 
non-Gaussian contributions are always negligible for TT 
and EE; for BB they are significant if the number of 
bands is small and Imax is sufficiently large, but can be 
hidden if the fit is performed using many bands. This 
is consistent with the discussion in non-Gaussianity 
is hidden when considering narrow I bands, but appears 
as extra variance when estimated BB power is averaged 
over a wide range in I. 



In Figure lTBl we have computed A(x^) using Eq. ljA3|l . 
and AVar(x^) using the approximation of Eq. (|A4() . To 
check this approximation, and the approximation that 
Afij can be computed to lowest order in Cf"^, we have 
also computed (x^), Var(x^) using Monte Carlo simula- 
tions of the lensed CMB, and find excellent agreement 
throughout Figure [TCI 



APPENDIX B: FIDUCIAL MODEL 
DEPENDENCE 

Throughout this paper, we have presented results for 
the fiducial model of Eq. (|23|l which has a low ioniza- 
tion optical depth and correspondingly a low erg = 0.73. 
For small deviations around the fiducial model, we have 
found that the shape of the principal components (Fig.|21) 
is unchanged, but the translation between the noise 
level and uncertainties a{Qi) in the lensing observables 
(Fig. EJ) can be affected. Denoting the uncertainty at 
noise level Ap by (j{Qi; Ap), we find the following rough 
scaling, expected from signal-to-noisc considerations as- 
suming that the unlensed CMB is fixed: 



a(ei;Ap) 
(T(e2;Ap) 



(j<i>4> 



Cf 



Cfid 



(Tfid(6i; Ap) 



(Bl) 



-l,fid, 



02; 



Ik2 



-1/2 ■ 

Ap 



/A'2.fid , 



Here, 
from 



Kl 



114, Ik2 = 440 are the median multipoles 
The scaling for Qi follows from considering 



the unlensed CMB as a fixed noise source whereas for 02 
it follows from direct signal-to-noise scaling. 

The optimal noise level from Eq. 147(1 for measuring 

lensing i3-modes scales roughly as {Cf^_^/Cf^^ fid)^^'^ ^^o'^ 
the same reason. 
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